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Abstract 

We investigate the behavior of an electric potential profile inside a mesoscopic 
region attached to a pair of superconducting leads. It turns out that the I — V 
characteristic curves are strongly modified by this profile. In addition, the 
electronic population in the mesoscopic region is affected by the profile behavior. 
We derive the single particle current and the mesoscopic electronic population 
within the non-equilibrium Keldysh Green functions technique. The Keldysh 
technique results are further converted in a self consistent field (SCF) problem 
by introducing potential profile modifications. Evaluation of I— V characteristics 
are presented for several values of the model parameters and comparison with 
current experimental results are discussed. 
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1. Introduction 

From the demostration of a superconductor-normal-superconductor (S-N- 
S) transistor [l|, the study of the nonequilibrium transport through supercon- 
ducting systems has been of much interest [1, H, 0, S Q • Another interesting 
problem in mesoscopic physics is transport through a superconductor/quantum 

dots/superconductor system @, H S El El, 0, El HQ- in this paper, we 

study the effect of an electrostatic potential profile on the electric transport 
across a single quantum dot with a spin degenerated level. Such a dot is cou- 
pled to a pair of biased superconductors contacts or leads (source and drain). 
By applying a source voltage Vs and a drain voltage Vd an electric current 
can flow between the leads and across the quantum dot which sets a typical 
non equilibrium situation. Besides the applied drain voltages Vd and source 
voltage Vs the system is further manipulated by a gate voltage Vg which, in 
principle, couples directly to the quantum dot. It turns out that Vd, Vs and Vg 
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induce an effective electrostatic profile potential inside the mesoscopic region 
in such a way that electronic population and electric current become tied to a 
self consistent problem. It is quite clear that such potential profile modifies the 
quantum dot level structure in a self consistent fashion. Such situation can be 
highly complicated since it mixes non equilibrium statistical mechanics with a 
classic electrostatic framework. Here, we adopt an approach which relates the 
self consistent electrostatic profile to the electronic population of the quantum 
dot and to the electric current [l6, 17 1. The self consistency and any other model 
calculations are fully performed within the non equilibrium Keldysh technique 

3 El- 
la section [5] we find the expression for the current and the electronic popu- 
lation for a mesoscopic region. In addition, we show calculations which lead to 
a self consistent field (SCF) problem between the dot electronic population and 
the electric current between the superconducting leads. The self consistency 
takes into account electric potential profiles inside the mesoscopic region as in- 
duced by the drain and source bias and by the gate voltage [17|. Moreover we 
present the results about the effect of the potential profile on the 1—V charac- 
teristic curves and on the electronic population inside the mesoscopic region. 
Finally, in section U we discuss our main conclusions. 



2. Calculation 

In this section we present the model and calculations which lead to the 
current and to the population number in the mesoscopic region. 

We consider a spin degenerated single orbital in a quantum dot connected 
to superconductors leads. The hamiltonian which describes this system is a 
generalized Fricdcl- Anderson model 20] . It reads 

H = H S + H D + H T , (1) 

where H$, Hjj and Ht stand for the superconducting leads, the dot and the 
tunneling term, respectively. H$ = -Hl + Hr where Hl and Hr are the left and 
right lead hamiltonians, respectively. They are given, within the BCS model 
El, by 

H s = V* f - H° r * t (2) 

rjka 

with 

ttO _ | £ vk ^r/k \ 

where is the superconductor gap, of the lead r\ 
the Nambu spinors 
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Hd is the hamiltonian for the single-level quantum dot of energy Ed'. 



(5) 



with 



R QD 



E d + U d n_ CT 

-E d - \J d n a 



The tunneling hamiltonian Ht is given by 



r/ka 



(6) 



(7) 



with 



(8) 



Ht connects the dot to the biased superconducting leads and it allows the elec- 
tric charge flow. VLg is the hybridization matrix element between a conduction 
electron of energy in the rj = L, R superconductor lead and a localized 
electron on the dot with energy <jy a and0<j are the dot spinors 



(4 d -c 



da 
4- 



0) 



Here cl 1 - fc^g^ denotes the creation (annihilation) operator for a conduction 

electron with the wave vector fc, spin a in the rj = L,R superconductor lead. 
d\ (d a ) is the creation (annihilation) operator for an electron on the dot. 
The flow of electric charge from the terminal rj is given by 



l v (i) = (-e) 



d(JV„(t)> 



dt 



[H T (t),N v (t)]) 



(10) 



where — e is the electron charge. (• • •) is the thermodynamical average over the 
biased L and R leads at the temperature T. Equation (fTU)) can be expressed in 
terms of the Keldysh Green function 



as 



F„jE« r (*.0 = -i<TcO njeir (t)dt(*')> 
2e. 



Ir, (*) 



(ii) 



(12) 



where (t,t') is a lesser Keldysh Green function. For the purpose of the 
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single particle current evaluation the coupling 



can be replaced by an 



average at the Fermi surfaces of the leads L and R. Hereafter, for simplicity, 
we replace V g by as we already do it in eqn ([1 
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The first evaluation step of (t, f) expresses it in terms of dot Keldysh 
Green functions. Then, we set an equation of motion for the Keldysh Green 
function F nKa (t,f) 

4 - *W (*. = ~^ v %a (*, *) + V V G„ (t, t>) , (13) 



where 



^ (*.*') = -i (TcaJ - 4(f)), (14) 

G CT (t,f) = -i(T c tk(i)4(*')>- (15) 
Similarly, .T 7 g CT (t, f ) satisfies the equation of motion 

(i| + Se) ^fc (*» = A ^ (*« ~ (*» - ( 16 ) 

where 

a CT (t,f) = -i(T c dL CT (t) 4 (*'))■ (i7) 

The eqns (fT3|) and (|16|) can be written as follows: 



[ -A, if + ^ J I ^ (t, f ) J " V- { g a (tj J ■ (18) 
This equation can be written as an integral along the Keldysh contour Ck 



; " i CK V W (*•*") (*•*") / 

W2'^$ ), (19) 



The 2x2 matrix in the right hand side of eqn (|19p is an unperturbed Keldysh 
Green function where 

g,fc (*'*') = -K^V-W^^*'))^ (20) 

(*.0 - - i ( Tc0 I.-s,- ( * )o «».-«.- (0 )o' (21) 



- -K Tc <-*v^H^) - (23) 

The subindex o indicates that evaluations are performed with V n = 0. 

The contribution (t, f ) S p to the single particle (SP) current is given by 
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F^Msp = v n / 

J — i. 



At" x 



;« ct (t, i") G< (f, + g<^ (t, f") G« (*", f )] (24) 



where we used eqn (fT9|) and Langreth rules 22 |- The superscripts "^^W'W 
correspond to lesser, greater, retarded and advanced Keldysh Green functions, 
respectively. Therefore, the single particle current I,, (t) SP can be written as 



I, (t) 



SP 



dt' 



E« (t, t>) G< (*', t) + E< CT (t, f ) G« (f , t)l (25) 



E$ (t,f) - Efeg^l (M') and E< (i,f) = V, 2 £gg< ^ (M') are self ener- 
gies which are evaluated for isolated superconductors leads L y R. They depend 
on t and t' through t — t' and are independent of a. Their Fourier transforms 
are given by 



= -r„ 



Mr, 



A, 



C( A r,, w ~ Mr,) + iC( w - Mr,, A^) 



where 



2T v C(lj - A,,)f(w - /i r) ) 



Vw 2 — cj' ; 



(26) 
(27) 

(28) 



r j; = 7riV^ (0) I^ 2 are the coupling constants between the leads and the quantum 
dot in the wide band limit. N n (0) is the density of states at the r/ Fermi level 
and f (lj) is the Fermi-Dirac distribution function. 
Equation (|25|) becomes 



Ir, (*) 



SP 



h 



du 



^e-K— ')« x 
2tt 



E« ( W ) ^ G< ( Wl «') + E< (c) £ ( w , «') 



(29) 



Dot Keldysh Green's functions G^ (lj,lj') and (lj,lj') are straightfor- 
ward evaluated. It turns out that they are a independent and frequency diagonal 
in the stationary limit 



G<(lj,lj') = 

G^(LJ,Lj') EE 



2ir6 (lj-lj')G < (lj) 
2it5 (lj - lj') G (a) (lj) . 



(30) 
(31) 



Equation (|29)) becomes 
4e 



I 



dw 



EWHG<M + S<HG< a 'H] (32) 
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The final expression for the single particle current Isp = (Ir,sp ~ ^l,sp) /2 
is given by 

lsp = -r- do; — — — — - p w f w - eV - f (w) , (33 

ft J-oo T L (w - eVj + (uj) 



In eqn (|33[) we performed a trivial shift of the dot energy level and insert 
the electric potencial V through eV — [J,l — fiR- The extra 2ir factor arises from 
the dot Keldysh Green functions. r r) (w) and p (u>) are given by 

r,» = r„c(w,A„) (34) 

PH = -1SKJW(«) = T -M^ 

{u-E d f+T*{u) 

r(w) = r L ( W - e y) + t r («) (35) 

Here p (w) is the so-called quantum dot spectral function which is given in 
terms of the retarded (w) Keldysh Green function [3]. At steady state there 
is no net flow into or out of the mesoscopic channel which yields a stationary 
particle number in it. The population number N, at the dot, is given by 

/OO 1 
^G<M, (36) 

which becomes a weighted average over the L and R contacts 



r (uj) r (u) 



(37) 



So far, we are not included the side effects of a potential profile inside the 
mesoscopic channel. Such potential is induced by the action of source, drain 
and gate applied voltages. Since the number of quantum levels in the channel is 
small the particle number variation is negligible. It amounts to neglect potential 
profile variation inside the channel. Then we can visualize the channel as a single 
point and an equivalent circuit framework is quite useful. In this framework we 
associate capacitances Cp> , C5 y Cq to the drain, source and gate, respectively. 
Whenever drain, source and gate bias potentials Vn, Vs and Vg, respectively, 
are present it induces a shift U = — e (V c h — Vo) of the electrostatic energy inside 
the channel. V c h and Vo are channel electrostatic potentials after and before we 
apply the source and drain biases, respectively. The electronic population before 
and after we apply the biases mentioned above are given by 

-eiV = C D V + C S V + C G V (38) 
-eN = C D (V ch -V D ) + C s (V ch -Vs) + C G (V ch -V G ), (39) 

respectively. It leads us to 

- eAiV = —e(N — N ) - C E (V ch - V ) - C D V D - C S V S - C G V G (40) 
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where Ce = Cp> + Cs + Cq. The electrostatic potential shift U inside the 
channel becomes 

e 2 

^ = C/£ + — AiV (41) 

where 

C E Uc = C D (-eVo) + C s [-eV s ) + C G (-eV G ) (42) 

The first term yields linear contributions to the potential profile while the second 
one introduces a direct dependence on the electronic population N. Uq = e 2 /Ce 
is the dot charging energy. Ce is an effective dot capacitance which depends 
on drain C'e>, source Cs and gate Cq capacitances within an equivalent circuit 
framework. 

The potential profile U shifts the dot quantum levels such that Isp and N 
are found from a system of self consistent equations. 

t 8ne f°° a r L (uj - eV)r R (uj) 

1 sp = —r- I dw — — tTTXt 777\ PK U ^ U ) f ( w ~eV)-i (w) , (43 

h J-oo T L (w - eV) + T R (w) 



j-oo r L (w - e^) + r fl (cj) 

Equation (1441) determines iV in a self consistent fashion which inmediately 
yields the single particle electric current Isp by carrying out the integration in 
eqn (j43]) . 

We will consider a situation where the lead couplings are not extremely small 
and the dot capacitance is reasonably large. It will smear out the Coulomb 
blockade effect and the double occupancy of the resonant level will be very 
unlikely. 



3. Results and Discussion 

In the dashed curve in Figure Q] we show zero temperature I — V characteris- 
tics, calculated without the self consistent field (non-SFC) method for values of 
gate voltage V g > 0. In the same figure the solid curve shows zero temperature 
I — V characteristics, calculated with the self consistent field (SCF) method for 
values of gate voltage V g > . As we can see in the dashed curve, the current 
is nonzero for positive values of the drain voltage, while for negative values of 
the drain voltage the current vanishes out. For the solid curve the current can 
have nonzero values. 

In Figure [5] we show zero temperature I — V characteristics for gate voltage 
values V g < which are calculated without the self consistent (non-SCF) method 
(dashed curve) and with the self consistent method (SCF) (solid curve). In the 
first case (non-SCF) the current is zero for positive values of the drain voltage 
and nonzero for negative values of the drain voltage, while n the second case 
(SCF) we can observe a symmetric I — V characteristic. 
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Figure 1: Zero temperature I-V characteristics calculated without the self consistent field (non- 
SCF) method (dashed curve) and calculated using the self consistent field (SCF) method (solid 



curve), with = 
0.008 meV and A 



0.2 meV, V a 
= 0.2 meV. 



1 meV, U = 0.0025 meV, C D /C E = 0.5, T L = F R 
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Figure 2: Zero temperature I-V characteristics calculated without the self consistent field (non- 
SCF) method (dashed curve) and calculated using the self consistent field (SCF) method (solid 
curve), with E d = 0.2 meV, V g = -1 meV, U = 0.0025 mcV, C D /C E = 0.5, T L = T R = 
0.008 meV and A = 0.2 meV. 



9 



On the other hand, when we use the self consistent field (SCF) method, the 
single particle current reaches the maximum value for higher drain voltage as 
compared to the non-SCF method. It means that the presence of the potential 
profile U inhibits the electron flow. Furthermore, it is noted that the I — V 
characteristics, as calculated with the self consistent field (SCF) method, agree 
with experimental data reported in the literature 0. In adittion, it proves 
that the single particle electric current can have nonzero values for positive and 
negative values of the drain voltage. 

4. Conclusions 

In conclusion, we have studied the effect of a self consistent potential profile 
on single particle electric current across a mesoscopic system attached to su- 
perconductor leads. Such system describes a spin degenerated single quantum 
where Coulomb blockade is neglected, in the regime A Tl^r. We derived 
an exact single particle electric current by means of the many body Keldysh 
technique. Zero temperature I — V characteristics agree with the experimen- 
tal results. Furthermore, we showed there are symmetric I — V characteristic, 
within the self consistent method which address an interplay among potential 
profile, single particle electric current and electronic population at the meso- 
scopic region. 

References 

[1] J. J. A. Baselmans, A. F. Morpurgo, B.J. Wees and T. M. Klapwijk, Nature 
397, 43 (1999). 

[2] A. F. Volkov, Phys. Rev. Lett. 74, 4730 (1995). 

[3] F. K. Wilhelm, G. Schon, A. D. Zaikin, Phys. Rev. Lett. 81, 1682 (1998). 
[4] S. K. Yip, Phys. Rev. Lett. 58, 5803(1998). 

[5] F. Giazotto, T. T. Heikkila, F. Taddei, R. Fazio, J. P. Pekola and F. Bel- 
tram, Phys. Rev. Lett. 92, 137001 (2004). 

[6] M. G. Pala, M. Governale and J. Koing, New J. Phys. 9, 278 (2007). 

[7] , S. De Franceschi, L. Kouwenhoven, C. Schonenberger and W. Wensdorfer, 
Nature Nanotechnology 5, 703 (2010). 

[8] R. Fazio and R. Raimondi, Phys. Rev. Lett. 80, 2913 (1999). 

[9] R. Fazio and R. Raimondi, Phys. Rev. Lett. 82, 4950 (1999). 
[10] K. Kang, Phys. Rev. B57, 11891 (1998). 
[11] P. Schwab and R. Raimondi, Phys. Rev. B59, 1637 (1999). 



10 



[12] A. A. Clerk, V. Ambegaokar, S. Hcrshficld, Phys. Rev. B61, 3555 (2000). 

[13] S. Shapira, E. H. Linfield, C. J. Lambert, R. Seviour, A. F. Volkov, A. V. 
Zaitsev, Phys. Rev. Lett. 84, 159 (2000). 

[14] Y. A. Levy, J. C. Cuevas, A. Lopez-Davalos, A. Martm-Rodero, Phys. Rev. 
B55, R6137 (1997). 

[15] J. C. Cuevas, Y. A. Levy, A. Lopez-Davalos, A. Martm-Rodero, Phys. Rev. 
B63, 094515 (2001). 

[16] G. C. Liang, A. W. Ghosh, M. Paulsson, S. Datta, Phys. Rev. B69, 115302 
(2004). 

[17] S. Datta, Electronic Transport in Mesoscopic Systems (Cambridge Univer- 
sity Press, UK, 1997). 

[18] L. Keldysh, Sov. Phys. JETP 20, 1018 (1965). 

[19] H. Haug and A. P. Jauho, Quantum Kinetics in Transport and Optics of 
Semiconductors (Springer, Berlin, 1996). 

[20] P. W. Anderson, Phys. Rev. 124, 41 (1961). 

[21] J. Bardcen, L. N. Cooper and J. R. Schricffcr, Phys. Rev. 108, 1175 (1957). 
[22] D. C. Langreth, Phys. Rev. 150, 516 (1966). 

[23] D. C. Ralph, C. T. Black and M. Tinkham, Phys. Rev. Lett. 74, 3241 
(1995). 



11 



